A Multiparameter, Numerical Stability Analysis of a Standing Cantilever Conveying Fluid

In this paper, we numerically examine the stability of a standing cantilever conveying fluid in a multiparameter space. Based on nonlinear beam theory, our mathematical model turns out to be replete with exciting behavior, some of which was totally unexpected and novel, and some of which confirm our intuition as well as the work of others. The numerical bifurcation results obtained from applying the Library of Continuation Algorithms (LOCA) reveal a plethora of one, two, and higher codimension bifurcations. For a vertical or standing cantilever beam, bifurcations to buckled solutions (via symmetry breaking) and oscillating solutions are detected as a function of gravity and the fluid-structure interaction. The unfolding of these results as a function of the orientation of the beam compared to gravity is also revealed

Scaling behavior in economics: I. Empirical results for company growth

We address the question of the growth of firm size. To this end, we analyze the Compustat data base comprising all publicly-traded United States manufacturing firms within the years 1974-1993. We find that the distribution of firm sizes remains stable for the 20 years we study, i.e., the mean value and standard deviation remain approximately constant. We study the distribution of sizes of the ``new'' companies in each year and find it to be well approximated by a log-normal. We find (i) the distribution of the logarithm of the growth rates, for a fixed growth period of one year, and for companies with approximately the same size $S$ displays an exponential form, and (ii) the fluctuations in the growth rates -- measured by the width of this distribution $\sigma_1$ -- scale as a power law with $S$, $\sigma_1\sim S^{-\beta}$. We find that the exponent $\beta$ takes the same value, within the error bars, for several measures of the size of a company. In particular, we obtain: $\beta=0.20\pm0.03$ for sales, $\beta=0.18\pm0.03$ for number of employees, $\beta=0.18\pm0.03$ for assets, $\beta=0.18\pm0.03$ for cost of goods sold, and $\beta=0.20\pm0.03$ for property, plant, & equipment.Comment: 16 pages LateX, RevTeX 3, 10 figures, to appear J. Phys. I France (April 1997

Scaling behavior in economics: II. Modeling of company growth

In the preceding paper we presented empirical results describing the growth of publicly-traded United States manufacturing firms within the years 1974--1993. Our results suggest that the data can be described by a scaling approach. Here, we propose models that may lead to some insight into these phenomena. First, we study a model in which the growth rate of a company is affected by a tendency to retain an ``optimal'' size. That model leads to an exponential distribution of the logarithm of the growth rate in agreement with the empirical results. Then, we study a hierarchical tree-like model of a company that enables us to relate the two parameters of the model to the exponent $\beta$, which describes the dependence of the standard deviation of the distribution of growth rates on size. We find that $\beta = -\ln \Pi / \ln z$, where $z$ defines the mean branching ratio of the hierarchical tree and $\Pi$ is the probability that the lower levels follow the policy of higher levels in the hierarchy. We also study the distribution of growth rates of this hierarchical model. We find that the distribution is consistent with the exponential form found empirically.Comment: 19 pages LateX, RevTeX 3, 6 figures, to appear J. Phys. I France (April 1997

Power Law Scaling for a System of Interacting Units with Complex Internal Structure

We study the dynamics of a system composed of interacting units each with a complex internal structure comprising many subunits. We consider the case in which each subunit grows in a multiplicative manner. We propose a model for such systems in which the interaction among the units is treated in a mean field approximation and the interaction among subunits is nonlinear. To test the model, we identify a large data base spanning 20 years, and find that the model correctly predicts a variety of empirical results.Comment: 4 pages with 4 postscript figures (uses Revtex 3.1, Latex2e, multicol.sty, epsf.sty and rotate.sty). Submitted to PR

A Model of Vertical Oligopolistic Competition

This paper develops a model of successive oligopolies with endogenous market entry, allowing for varying degrees of product differentiation and entry costs in both markets. Our analysis shows that the downstream conditions dominate the overall profitability of the two-tier structure while the upstream conditions mainly affect the distribution of profits. We compare the welfare effects of upstream versus downstream deregulation policies and show that the impact of deregulation may be overvalued when ignoring feedback effects from the other market. Furthermore, we analyze how different forms of vertical restraints influence the endogenous market structure and show when they are welfare enhancing

HOMMEXX 1.0: a performance-portable atmospheric dynamical core for the Energy Exascale Earth System Model

We present an architecture-portable and performant implementation of the atmospheric dynamical core (High-Order Methods Modeling Environment, HOMME) of the Energy Exascale Earth System Model (E3SM). The original Fortran implementation is highly performant and scalable on conventional architectures using the Message Passing Interface (MPI) and Open MultiProcessor (OpenMP) programming models. We rewrite the model in C++ and use the Kokkos library to express on-node parallelism in a largely architecture-independent implementation. Kokkos provides an abstraction of a compute node or device, layout-polymorphic multidimensional arrays, and parallel execution constructs. The new implementation achieves the same or better performance on conventional multicore computers and is portable to GPUs. We present performance data for the original and new implementations on multiple platforms, on up to 5400 compute nodes, and study several aspects of the single- and multi-node performance characteristics of the new implementation on conventional CPU (e.g., Intel Xeon), many core CPU (e.g., Intel Xeon Phi Knights Landing), and Nvidia V100 GPU.</p

A mouse chromosome 4 balancer ENU-mutagenesis screen isolates eleven lethal lines.

BACKGROUND: ENU-mutagenesis is a powerful technique to identify genes regulating mammalian development. To functionally annotate the distal region of mouse chromosome 4, we performed an ENU-mutagenesis screen using a balancer chromosome targeted to this region of the genome. RESULTS: We isolated 11 lethal lines that map to the region of chromosome 4 between D4Mit117 and D4Mit281. These lines form 10 complementation groups. The majority of lines die during embryonic development between E5.5 and E12.5 and display defects in gastrulation, cardiac development, and craniofacial development. One line displayed postnatal lethality and neurological defects, including ataxia and seizures. CONCLUSION: These eleven mutants allow us to query gene function within the distal region of mouse chromosome 4 and demonstrate that new mouse models of mammalian developmental defects can easily and quickly be generated and mapped with the use of ENU-mutagenesis in combination with balancer chromosomes. The low number of mutations isolated in this screen compared with other balancer chromosome screens indicates that the functions of genes in different regions of the genome vary widely

Capillary Condensation and Interface Structure of a Model Colloid-Polymer Mixture in a Porous Medium

We consider the Asakura-Oosawa model of hard sphere colloids and ideal polymers in contact with a porous matrix modeled by immobilized configurations of hard spheres. For this ternary mixture a fundamental measure density functional theory is employed, where the matrix particles are quenched and the colloids and polymers are annealed, i.e. allowed to equilibrate. We study capillary condensation of the mixture in a tiny sample of matrix as well as demixing and the fluid-fluid interface inside a bulk matrix. Density profiles normal to the interface and surface tensions are calculated and compared to the case without matrix. Two kinds of matrices are considered: (i) colloid-sized matrix particles at low packing fractions and (ii) large matrix particles at high packing fractions. These two cases show fundamentally different behavior and should both be experimentally realizable. Furthermore, we argue that capillary condensation of a colloidal suspension could be experimentally accessible. We find that in case (ii), even at high packing fractions, the main effect of the matrix is to exclude volume and, to high accuracy, the results can be mapped onto those of the same system without matrix via a simple rescaling.Comment: 12 pages, 9 figures, submitted to PR

Albany: Using Component-based Design to Develop a Flexible, Generic Multiphysics Analysis Code

Abstract: Albany is a multiphysics code constructed by assembling a set of reusable, general components. It is an implicit, unstructured grid finite element code that hosts a set of advanced features that are readily combined within a single analysis run. Albany uses template-based generic programming methods to provide extensibility and flexibility; it employs a generic residual evaluation interface to support the easy addition and modification of physics. This interface is coupled to powerful automatic differentiation utilities that are used to implement efficient nonlinear solvers and preconditioners, and also to enable sensitivity analysis and embedded uncertainty quantification capabilities as part of the forward solve. The flexible application programming interfaces in Albany couple to two different adaptive mesh libraries; it internally employs generic integration machinery that supports tetrahedral, hexahedral, and hybrid meshes of user specified order. We present the overall design of Albany, and focus on the specifics of the integration of many of its advanced features. As Albany and the components that form it are openly available on the internet, it is our goal that the reader might find some of the design concepts useful in their own work. Albany results in a code that enables the rapid development of parallel, numerically efficient multiphysics software tools. In discussing the features and details of the integration of many of the components involved, we show the reader the wide variety of solution components that are available and what is possible when they are combined within a simulation capability. Key Words: partial differential equations, finite element analysis, template-based generic programmin